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Abstract 

Using the O(N) model as an example, we investigate the self-interaction effects of inflaton on 
the dynamics of the primordial perturbations. When taking interactions into account, it is essen- 
tial to employ a self-consistent off-equilibrium formalism to study the evolution of the inflationary 
background field and its fluctuations with the back-reaction effects. Within the Hartree factoriza- 
tion scheme, we show that the 0(N) model has at least two observable remains left behind the 
off-equilibrium processes: the running spectral index of primordial density perturbations and the 
correlations between perturbation modes in phase space. We find that the running of the spectral 
index is fully determined by the rate of the energy transfer from the inflationary background field 
to its fluctuations via particle creation processes as well as the dynamics of the background field 
itself. Furthermore, the amplitude of the field fluctuations turns out to be scale-dependent due 
to the off-equilibrium evolution. As a consequence, the scale-dependence of fluctuations yields a 
correlation between the phase space modes of energy density perturbations, while the one-point 
function of the fluctuations in each Hartree mode is still Gaussian. More importantly, the mode- 
mode correlation of the primordial perturbations depends upon the dynamics of the self-interaction 
as well as the initial conditions of the inflation. Hence, we propose that the running spectral index 
and the correlation between phase-space modes would be two observable fossils to probe the epoch 
of inflation, even beyond. 
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I. INTRODUCTION 



In the inflationary scenario, the primordial perturbations of the universe originate from 
vacuum fluctuations of the scalar field(s), the inflaton <f), driving the inflation. If the dynam- 
ics of the fluctuations is approximated by single massless free field during the inflationary 
epoch, the power spectrum of curvature perturbations for a Fourier mode k can be obtained 

as 



Vn{k) = 




where the inflationary scale factor a = exp (Ht), and the Hubble parameter H = ^j8ttGV/3 

are determined by the potential of the inflaton field V(<p). The first factor H/(j) on the 
right hand side of Eq. (Q) comes from the evolution of the background inflaton field 0, i.e. 
the expectation value of the quantum scalar field. Meanwhile, the second factor H/2tt is 
specified by the variance of the classical <fi field fluctuations, 5(pk, at a few Hubble times after 
horizon-crossing. With slow-roll conditions, the derivative of the spectral index, dn(k) /dink, 
turns out to be very small, or even negligible such that n(k) = (d\n.Vn{k) / dhik) +1 — 1 is 
roughly fc-independent 

However the fc-dependence, or the running of the spectral index, revealed by the re- 
cently released Wilkinson Microwave Anisotropy Probe (WMAP) data can be as large 
as dn(k)/d\nk = (d/d\nk) 2 \nV n (k) ~ -0.055 to - 0.077 [2]. Hence, the ori- 
gin of the primordial perturbations cannot be solely accommodated with the quan- 
tum fluctuations of single free field. In a slow-roll inflation, one has d/dlnk = 
(1/ ' H)d/dt = —( y l/SiTG){V'((j))/V((f)))(d/d(f)), and therefore the derivative of the power spec- 
trum (d/d In k) n In V-n(k) will no longer be negligible for the second order (n = 2) as long 
as the third or higher derivatives of the inflaton potential becomes substantial. To fit in 
with the running of the spectral index, models beyond single scalar field with quadratic 
potential V(4>) have been proposed accordingly In this context, the running index of the 
perturbation power spectrum is considered as an essential feature due to the interactions or 
the self- interactions of the inflaton(s). 

It should be pointed out here that in deriving the derivative of the power spectrum 
the tree level of the effective potential V(<f)) is implicitly involved. However, when quantum 
fluctuations arising from the interactions of the inflaton become important, this approach 
may be problematic. In particular, the one-loop effective potential will turn complex within 
the region where the background field (ft is constrained by V"((f>) < 4]. The imaginary 
part of the effective potential would inevitably lead to a dynamically unstable state p. 
This so-called "spinodal instability" will allow long wavelength fluctuations to grow non- 
perturbatively p. As such, the term H/2% in Eq. (1) is no longer valid when the inflaton 
suffers from spinodal instabilities. Therefore, the primordial perturbations must be dealt 
consistently with a different method to account for the amplification of vacuum fluctuations 
in the presence of spinodal instabilities. 

This motivates us to study the self-interaction effects of the inflaton on the dynamics of 
primordial perturbations within a context of the self-consistent off-equilibrium formalism. 
We will use the O(N) model of inflation as an example. When N = 1, it reduces to the pop- 
ular X(j) 4 inflationary model which is disfavored by the WMAP. However, the O(N) model 
with spontaneously broken symmetry is still a viable model in quantum field theory and 
off-equilibrium statistical mechanics pj. It has been extensively used in modeling the quan- 
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turn off-equilibrium processes in the early universe, as well as the chiral phase transition in 
relativistic heavy ion collision, etc |8j. In particular, the dynamics of quantum fluctuations 
of the OJN) model in the large N limit has been developed with the Hartree-type factoriza- 
tion which can be generalized for the finite N situation. In this paper, we will focus 
on searching for the possible observable imprints caused by the off-equilibrium evolution of 
inflationary primordial perturbations. 

The paper is organized as follows. In Section 2, we briefly describe the O(N) model of 
inflation, and introduce the self-consistent, nonperturbative and renormalized solutions to 
the nonlinear evolution of the background inflaton field as well as its quantum fluctuations 
within the Hartree approximation. Section 3 analyzes the dynamics of quantum fluctua- 
tions undergoing a quantum-to-classical transition, and the statistical properties of the cor- 
responding classical primordial perturbations. The spectral index along with the running of 
the perturbation power spectrum will be calculated in Section 4. Section 5 then discusses 
the correlation of primordial perturbations in phase space and its detection. Finally, we 
summarize our findings and give conclusions in Section 6. 



II. THE 0(N) INFLATION 



A. The model and the Hartree factorization 

Consider the dynamics of an inflation driven by a field <fr of the O(N) vector model with 
spontaneous symmetry breaking. The action is defined by 



-9 



S = J d 4 xC = J d A x 
where V(*& ■ is a self-interaction potential given by 

V(* ■ = — I # • * 



_A_ 



2Nm 2 \ 2 



A 



(2) 



(3) 



As an inflaton, the N components of the field generally are represented as $ = (a, ff), where 
7? represent N — 1 scalar fields. The cosmic inflation is characterized by the state in which 
the component a has a spatially homogeneous expectation, i.e. a can be decomposed into a 
background plus fluctuations around the background as 

a(x,t) = ^/N(f>(t)+x(x,t) (4) 

with the expectation value (cr(x, £)) = \/N(p(t) and thus, (x(x, £)) = 0. During the infla- 
tionary epoch, the background space-time can be described by a spatially flat Friedmann- 
Robert son- Walker metric, 



ds 2 = g^dx^dx 11 = dt 2 — a 2 {t)5ijdx l dx^ , 

where the scale factor a = exp(Ht) with the expansion rate H = 
the mean of energy density of the inflaton field as 



(5) 



H 2 (t) 



8nN 
3M 2 X 



1 • 



A 



4>\t) + - <f>\t) 



2m 



2\ 2 ' 



8irGp/3 determined by 
(6) 
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where Mp\ is the Planck mass. We have assumed that the effects of quantum fluctuations 
on the dynamics of the Hubble parameter can be ignored as we will justify this assumption 
later. 

During the inflation, the expectation value of the scalar field <f>(t) undergoes an off- 
equilibrium evolution from the initial state <p(t) ~ where V"(4>) < 0. Accordingly, the 
mass square of the long wavelength fluctuation modes will be negative, which leads to 
the nonperturbative growth of fluctuations. Therefore, a nonperturbative framework is 
necessary for taking account of the growth of fluctuations, especially in computing the power 
spectrum of perturbations later. We will employ the method of the Hartree factorization, 
which approximates the potential V(4>) with an effective quadratic potential while keeping 
N finite 0, ■ The Hartree- factorized Lagrangian is 



C(t) = / d 3 xa 3 (t) 



+ 2 0m*) 



\Ml{t) X * - l -Ml{t)^ 



xVit) 



where 



V'(t) = VN y>(t) + 3Hj>(t) + {M\ - A0 2 (t))0(t) 
The time dependent effective masses M x (t) and M n (t) are obtained as 



(7) 



(8) 



M 2 x {t) 



M 2 (t) 



-m 2 + y0 2 W + |^(x 2 )(t) 



~m 2 + \m + ±tf)(t) + ± (l + ^ ) (^ 2 )(t) 



(9) 



where (ip 2 ) is defined by (n 2 ) = (JV — 1)(V ;2 )- 

The Hartree approximation we adopt here is equivalent to the Hartree-Fock approxima- 
tion in the general CJT (Cornwall, Jackiw and Tomboulis) formalism 0. By the requirement 
of Hartree approximation, the one-loop corrections of the quantum field two-point Green's 
functions are cancelled by the introduced "mass counter term". For the case of finite N, 
however, the 1/N corrections in the Hartree factorization do not include the contributions 
from the collision effects of the same order of 1/iV HJ. 

Nevertheless, the factorization 
method provides a reliable resummation scheme that allows us to treat the growth of quan- 
tum fluctuations driven by spinodal instabilities self-consistently. 

The equation of motion of the background field <fi can be directly obtained from the term 
xV{t) i n t ne Lagrangian [Eq. (JJJ)]. By means of the tadpole condition (x(x, £)} = and 
Eq. (jSJ), we have 



<j>{t)+3H(f>{t)+ Ml(t)-\<p 2 (t) <p(t) = 0. 



(10) 



To see how the quantum effects influence the dynamics of the background field, it is critical 
to solve the background field 4>(t) self-consistently by including the fluctuations (x 2 )(t) and 
(^ 2 )(t) in the mass-squared term M 2 (t). 

To find the dynamics of the fluctuation fields, we decompose x(t) and n(t) in the Fourier 
basis. In the Heisenberg picture, one has 



xCM) 

7Ti(x,t) 



d 3 k 
8^ 
d 3 k 



x(M) : 

%i(k,t) 



d 3 k r 



b k f x , k (t) + &L£,_ k (*) 



„ikx 



d 3 k 



(11) 
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where a^, and aj k , b\. are the creation and annihilation operators which obey the commu- 
tation relations [ct^, a] k ,] = 5^8^ and [6k, b\,,] = 8 k k '. The equations of the mode functions 
/x,k(£) and fir^it) can be found from the Heisenberg field equations given by 



fx*® = 0, 
U*(t) = 0. 



(12) 



Finally, to close these equations self-consistently, the terms (x 2 )(t) an d (i/j 2 )(t) in the 
mass-squared [Eq. (JJJJ)] can be determined by the mode functions as 

d 3 k , 



87T 



"l/x.k(*)l' 

IA,k(t)P 



(13) 



Obviously, fluctuations (x 2 )(t) and (ip 2 )(t) are created from the self-consistent background 
state 4>(t) during the inflation. The nonlinearity, or the self-interaction, is in fact encoded 
in the self-consistent solutions to Eqs. ((TUJ) - ((TSJ) . 



B. Renormalization and non-equilibrium equations of motion 



Before proceeding further, one needs to understand the issue of renormalization regarding 
the divergences associated with the loop integrals. The divergences can be determined from 
the large loop momentum behavior of the mode functions, which can be found from the 
WKB-type solutions to Eq. (|T2~j) 10]. It turns out that the self-consistent loop integrals 
in Eq. (fT3*j) contain both quadratic and logarithmic divergences. To get the self-consistent 
renormalized equations, one has to subtract the divergences from both the bare mass and the 
coupling constant. However, for a weak coupling A < 10~ 14 in a typical inflation model, the 
logarithmic subtractions can be neglected. Thus, the mass renormalization can be simplified 
as 

1 A 2 " 



m. 



m 



A / 2 
+ 2V + N 



in 2 a 2 



with a negligible renormalization for the coupling constant, 
renormalized masses M Xt R and M % ,r which are given by 



(14) 

As a result, it leads to the 



A 



Ml R {t) 



A 



A 



-m, 



+ ~r(t) + —( x 2 }R(t) + - 



where 



(X 2 )R(t) 



2N 



A d 3 k 



'4 



(15) 



57T 



A d 3 k 



i/x,kwr- 
lum 2 - 



1_A^ 

]tt 2 a 2 
1 A 2 



Itt 2 a 2 



(16) 



Hereafter, we will drop the subscript R, as the integrals of mode functions are always in 
terms of the renormalized quantities. 
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III. CLASSICAL FLUCTUATIONS OF THE INFLATON 



A. Quantum decoherence 

The quantum-to-classical transition of the quantum fluctuations during the inflation can 
be investigated from the Schrodinger wave function approach jTij. [l3| . The Hamiltonian 
describing the evolution of quantum wave functions of the fluctuations can be obtained from 
Eq.©. In general, the developments of different fluctuation modes are separable due to the 
quadratic feature of the Hartree factorized Hamiltonian. Let us first consider a mode k of 
the x field. In the Schrodinger picture, the initial (time to) vacuum state is specified by 

bMs = 0, (17) 

for all k. At time t, the evolved Schrodinger state \t)s is given by 

U(t,t Q )b k U-\t,t o )\t) s = 0, (18) 

where U(t,t ) = exp —(i/h) TC(t)dt. In the coordinate representation, the conjugate mo- 
mentum can be written as II(k) = — ihd/dx(— k) , and thus, Eq.ffTSj) reduces to 



d X (-k) 



*[*(k),t]=0, (19) 



where 



^^nuW 11 - 2 ^-^' (20) 

F --M = ^5l/*, k (f)l 2 , (21) 

and ^[^(k),^ = (x(k)\t) s is the wave function for the field x(k). The solution to the 
equation can be obtained straightforwardly as 

*[*(k),t] = A/- k 1/2 («)exp (-^WMk)! 2 ) (22) 

= < 2(t)eX K-S2SF |1 " !2F '■ k(t,, )' 

where M^ 2 {t) is the normalization coefficient. The wave function of the quantum fluctua- 
tions is then the direct product of the wave functions for all k modes of the x field. Notice 
that for an initial vacuum state which is a pure state, it will remain the pure state under the 
unitary evolution. Therefore, the density matrix in the coordinate representation becomes 
the wave function times its complex conjugate given by 

np(x(k),x(k),t) = n*[x(k),*]**[x(k),*] 



tt Kf /,\ j X(k)x(~k) +x(k)x(-k) ' 
exp I ^-^yp [1 - z2F x>k (t)) 



(23) 



The density matrix for each mode k can be decomposed into the diagonal and off-diagonal 
elements as follows: 
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p[ X (k) + S(k),x(k)-S(k),t] 



(24) 



Afk(t) exp 



\x(k)\ 2 - ^F xM (t) fc(k)*(-k) + x(-k)*(k)] + |5(k)p 



One can immediately recognize that F x ^(t) is the phase of the off-diagonal elements. A 
quantum-to-classical transition is implied if the following condition is satisfied: 

F xM (t) > 1 . (25) 

Similarly for mode 7r, the condition of the quantum-to-classical transition is given by 



a 3 (t) d 
2 dl 



l/ 7r ,k(t)| 2 »l. 



(26) 



Hence, the quantum-to-classical transition is determined by the behavior of the time- 
dependence of the mode functions / x ,k(t) and f^^it). 

To solve Eq. ()12|). we consider the early stage of the inflation when fluctuations have 
not grown significantly yet in the presence of spinodal instabilities as well as the Hubble 
parameter remains constant. In this stage, the masses M x (t) and M n (t) can be approximated 
by —m in Eq. ()15|) . Thus, both mode equations in Eq. (jl2j) can be written as 



d 2 ^ H d k 2 
dt 2 dt a 2 



m 



h(t) = . (27) 
The general solution can be expressed as a combination of the Hankel function H u (k/aH) and 



the Neumann function N v {k/aH) with v = J (9/4) + (m/H) 2 . In particular, for fluctuations 
in superhorizon regime where k <C aH, the growing mode solution leads to 



fk ^ fo{k) e 



(v-3/2)Ht 



(28) 



where fo(k) depends on the initial conditions of the mode functions as we set to = 0. 
Apparently, F^t) 3> 1 is valid after several e-foldings of inflation when the modes become 
superhorizon-sized. 



B. Equations of classical fluctuations 

Using the density matrix, one can find the Wigner function from 



S(k) 



W( X (k) , 7r x (k) ,t) = JdU^j e-KW'OO p ^( k ) - x(k) + t 



P[|x(k)|] iV k exp 



where 



p[\m\] 



h 



%( k ) 



exp 



lx(k)| s 



K\fxAt)\ 2 \ 



(29) 
(30) 



Since the Wigner function in Eq.(J2HJ) is definite positive, it can be interpreted as a distribu- 
tion function in the phase space of stochastic field perturbations. 

For the perturbation modes just crossing out the horizon, the dynamics of the classical 
fluctuations can be described by the Wigner function in the limit of H — > while keeping 
^|/x>k( r )| 2 Axed. In fact, the term fi\f x ,u( T )\ 2 basically measures the variance of the field 
fluctuations x(k) and will grow nonperturbatively due to the spinodal instabilities. As a 
result, the Wigner function becomes 



W(x(k),7r x (k),t)=P[|x(k)|]5 



%(k) - F - k(t) 



l/x,k(*)l ; 



(31) 



The delta function on the above expression yields the equations of motion for the classical 
fluctuations at superhorizon scales, 



In terms of the classical fluctuation x(k, t), Eq. (|32jl can be rewritten as 



(32) 



(33) 



With the mode functions f x ,u.(t) derived self-consistently from Eqs. (11), (14), and (15), Eq. 
()33|) describes the dynamical evolution of the mean value of the stochastic fields x(k, t) in 
the superhorizon regime. Obviously, the stationary solution to the probability distribution 
function (PDF) of variables x(k, t) can be obtained from Eq. ([30)1 . 

In contrast to the usual second order evolution equation for the field x(k, t), Eq. (j22J) is 
a diffusion type equation which typically occurs during the overdamped process of harmonic 
oscillators. To accommodate the overdamped behavior of the mean field x(k, t) into the 
mode equations, the equations for the mode function Eq. (|T2~j) at superhorizon scales can be 
further approximated by 



3ff4 + -£r + M^(t) 



dt a 2 (t) 
As a consequence, Eq. (f3~3j) reduces to 



3H d * {Kt) + 



/*k(*) = o • 



(34) 



dt 



a 2 (t) 



+ M 2 (t) 



X(k,t) = 0. 



(35) 



As expected, these equations are the classical Euler-Lagrange field equations for the super- 
horizon modes of fluctuations x(k, t) in the overdamped regime. 

Following the same procedure, we can obtain a similar equation for the fluctuating field 
7r(k, t) as 

' ? ~ +A^(t)]7f i (k,t)=0. (36) 



3g^ (M) + 



dt 



a 2 (t) 



S 



C. The PDF of classical fluctuations 



The statistical properties of the classical fluctuation on mode k can be characterized by 
the PDF given by Eq. (J3Uj), 



P\k,(k) 



p[\xm 



i 



1/2 
I ex P 


[ ^ 2 (k) 


L l/ir,k(t)| 2 J 


v 1/2 

J exp 


X 2 (k) l 


L i/x, k (oi 2 J 



(37) 



Apparently, the PDFs of the classical perturbation modes 7Tj(k) and x(k, t) are Gaussian 
with variances \f-n-,k(t)\ 2 and |/ x ,k(i)| 2 respectively under the Hartree approximation. It 
implies that 

( X 2n+1 (k,t)) = 0; (^ +1 (k,t)) = 0, (38) 
( X 2 "(k,t)) = (2n-l)!!( X 2 (k,t))"; (^(k,*)) = (2n - l)!!<7r 2 (k, t)) n , 



where n = 1, 2, .... However, the variances generally are k-dependent. Such k-dependence of 
variances eventually leads to non-trivial correlations between fluctuation modes which will 
be scrutinized later. 

In summary, after a quantum-to-classical transition, the dynamics of the classical stochas- 
tic fields x(k, t), vfj(k, t) must be based on Eq. (J35j) and Eq. (J36)) with the evolution of the 
self-consistent mode functions. The corresponding PDF of the perturbations is governed 
by Eq. (j37J). In fact, by adding a noise term to the right hand side of both Eqs. (J55J) 
and (|36|). they become Langevin-like equations which is capable of describing the stochastic 
properties of perturbations. Thus, the noise correlations can be determined by the variances 
of the PDFs. For the case of single free inflaton, the dynamics of the super horizon- sized 
fluctuations is governed by a Langevin equation Il3j, which is obtained by coarse graining 
the degrees of freedom of all subhorizon modes [1^. According to the results we obtained 
here, the coarse grained Langevin equation approach is effective not only for the situation 
of a free scalar inflaton field, but also for models with self- interaction such as the O(N) 
inflation. 

Moreover, it should be pointed out that Eqs. (f3Ij) - (|35p and their counterparts for the 
fluctuating 7r(k, t)-field are not closed if considering only the classical fluctuations or the 
superhorizon modes. By virtue of M%(t) and M%(t), the subhorizon modes also contribute 
towards these diffusion type equations. This feature is anticipated because of the inevitable 
coupling between all wavelength modes under a substantial self-interaction. Thus, the renor- 
malized mass squared terms play critical roles as phenomenological potentials. Consequently, 
Eqs. (I34J) - (J33j) and their counterparts for the 7r(k, t) -field provide a possible scheme to 
study the interactions of inflaton by comparing the phenomenological potentials given by 
data fitting with the predicted M^(t) and M%(t). This characteristics is a main result of 
this paper. 

We now turn our attention to the cosmological implications from the off-equilibrium 
dynamics of the O(N) inflationary model in the next section. 
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IV. THE POWER SPECTRUM OF PRIMORDIAL PERTURBATIONS 



A. The mass density perturbations 

The power spectrum of primordial perturbations is described by V(k) = (\Su\ 2 ) where 
the mass density perturbations are determined by the gauge invariant quantity [ItJ 



5p 



P + P 



(39) 



k=aH 



The mass density fluctuation Sp originates from the field fluctuations 7Tj and x given by 



*! = (i-i 

N \ N 



with 1 



\{i> 2 ) + ^<W) 2 } - \m*m + ^ 2 «> 2 > + ^{xW) (40) 



+ ^ (i + ^ ) m 



2\2 



+ 



N 



^(x 2 ) + ^((Vx) 2 )-^ 2 (x 2 ) 



+ > 2 )^(i-i)(^)(x 2 }^<x 2 ) 2 



(x 2 )(t) 



d 3 k ■ 

\U(t)\ 2 



(2tt) 3 
d 3 k 
(2^)3 



\fxAt)\ 2 - 



1 A 4 

1 A 4 
37T 2 a 4 



+ 



+ 



^7r 2 a 2 
J_A^ 
37T 2 a 2 



M 2 (t) 



6 
R 

6 a 2 



(41) 



(W) 2 )(t) 
((Vx) 2 )(t) 



rf 3 k 



* 2 |/ir,k(*)| : 



1 A 4 A 2 



3-7T 2 a 2 



1 A 4 1A 2 / . R\ 
Mlit)--) 



Itt 2 a 2 8tt 2 V~ _7rw QJ 



(42) 



where the Ricci scalar R = 6(a 2 /a + a/a) = 12H 2 . Similar to the renormalization of Eq. 
(J16|) . we have ignored the logarithmic subtractions in the renormalization of above equations. 
The term summing up the energy density and the pressure, p+p, in Eq. can be obtained 
by 



P + P 
N 



1 

N 



<^ 2 > + -r< (W0 2 > 



1 

+ N 



(x 2 ) + -((Vx) 2 ) 



(43) 



Apparently, the first term 2 in Eq. (|4*3*|) comes from the background inflaton field. The 
other terms, however, are the contributions from the fluctuations. 



B. Numerical examples 

We present numerical examples in this section to demonstrate the self-interaction effect 
of the 4>{t) field upon the power spectrum. To begin with, we need to find the self-consistent 
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FIG. 1: The evolution of h = g 
for initial conditions -ff(O) = 2m, 0(0) 
case of iV = 4. 



2 m ' 2m 



(^), and (x 2 
0, 0(0) = 0, M 2 (0) = 



respectively vs. i (in units of m 1 ) 
M 2 (0) = m 2 with A = 10~ 14 in the 



solutions to the background field and the mode functions f Xt k(t) and f^^t) from Eqs. 
(fTU|) and (fT2j) . We assume 0(0) ~ and 0(0) = for the background field at the onset of the 
inflation. The initial (t = 0) conditions of the mode functions can be specified as effective 
free massive scalar fields in an expanding universe: 



/*,k(0) 

A, k (o) 



l 



2[F + M 2 (0)] 
1 

2[fc 2 + Af2(0)] 



, /x-k(0) 



, A,k(o) 



'A; 2 + M 2 (0) 



'fc 2 + M 2 (0) 



(44) 



where we have set a(0) = 1, and M 2 (0), M 2 (0) are the effective mass square of the field 
components x an d 71 , respectively. The values of M 2 (0) and M 2 (0) can be zero or a positive 
number, which depends on the details of the onset of the inflation. Although this uncertainty 
really is not decisive to the evolution of the mode functions / x ,k(*) and f nt k(t), it turns out 
to become crucial for correlating the fluctuations as we will see later. 

Figure 1 plots the expansion rate h(t) of the universe, the self-consistent solution to 
the background field (pit), and the fluctuations (ip 2 ){t), (x 2 )(£) where the parameters are 
taken to be A = 10~ 14 , N = A, and H{0) = 2m, M 2 (0) = M 2 (0) ~ m 2 > 0. The value 
for H(0) is determined from the energy density of the inflaton field given by Eq. @. 
During the inflation, the energy density of the inflaton is dominated by that of the potential 
energy which leads to h 2 (0) = (H/m) 2 w Nm 2 /M^X w 0(1) for N w 0(1), A w 10" 14 , 
and the inflaton mass m ~ 10 12 GeV. As a consequence, the Hubble parameter remains 
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FIG. 2: Mode function \f n ,]i(t)\ 2 vs. k (in units of m) at t = 40, 55, 70 (in units of m 1 ) under 
the same initial conditions and parameters as in Fig. 1 for the case of N = 4. 



approximately a constant h ~ 2 at t < 60m , but drops dramatically after t > 70m -1 . 
Apparently, t < 70m -1 specifies the inflationary epoch. 

It has been pointed out earlier that we treat the dynamics of the Hubble parameter 
H(t) classically without any quantum correction involved in Eq. @. It turns out that the 
background field does follow its classical equation of motion as the quantum corrections in 
Eq. (IK)jl give rise to a small Sp/p ~ 10 -4 in our case. Thus, when the semi-classical Einstein 
equations are employed for a truly self-consistent approach, it yields a similar solution of h 
as shown in [10]. 

Figure 1 also shows that <f>(t) ~ (ip 2 )(t) ~ (x 2 ){t) — during the first stage of the 
inflation t < 40m -1 in the spinodal regime where the initially positive M 2 and M 2 turn 
negative quickly. The spinodal instabilities become important eventually, and leads to the 
significant growth of both (ip 2 )(t) and (x 2 )(t) starting at t — 55m -1 . As expected, these 
properties are not sensitive to the choice of the initial effective masses M 2 (0) and M 2 (0). It 
can also be seen from the mode function |/,r,k(^)| 2 shown in Fig. 2. 

Evidently, the mode function undergoes a substantial increase in the period of 55m _1 < 
t < 70m -1 , but varies much slower after t = 70m -1 . This is owing to the predomination 
of the criterion M 2 < in the spinodal regions spawned during the slow rolling of the 
inflaton. On the other hand, the spinodal condition is only weakly satisfied over a span of 
70m -1 < t < 80m _1 where the background starts the rapid falling into the valley of the 
inflaton potential. Therefore, the process of spinodal instability terminates at time t e given 
by 

A , a,. A,n,., A/, 1 



Mi(t,) = -m' + -4,'(t, 



+ 



N 



(PKu 



o 



(45) 



It renders t e ~ (70 — 80)m 
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FIG. 3: The power spectrum of primordial density perturbations vs. In k (k in units of m ) with the 
same initial conditions and parameters as in Fig. 1 for the case of iV = 4, 20, and oo respectively. 



The sharp decline of h signifies the end of inflation. Therefore, the number of the total 
inflationary e-foldings is iV e « Ht = 2mt « 140. The perturbation just about to enter the 
present Hubble horizon should be transgressing the horizon in the inflationary era around 
t = t^o — 40m -1 , some 60 e- folds before the end of inflation. We find that most of energy 
density of the background (p state is transferred to the fluctuations (if) 2 ) via the production 
of 7i modes during the inflationary epoch. Since M 2 (t e ) ~ 0, these 7r modes are in fact 
the massless Goldstone modes of the broken symmetry phase. The Goldstone theorem is 
fulfilled dynamically. 

Using the numerical solutions, we calculate the power spectrum of primordial perturba- 
tions. The result is shown in Fig. 3, in which N is taken to be 4, 20, and oo. The N depen- 
dence of the power spectrum is trivial. We find that the power is generally ^-independent for 
modes k which cross the horizon when t < 55m -1 . While the energy transfer from the back- 
ground field to the fluctuations driven by the spinodal instability leads to a swift increase 
in (ip 2 ), the inflaton field at later times rolls down the potential hill which also results in 
an increase of 0. This in turn renders the gradual decrease in power for shorter wavelength 
modes crossing out of the Hubble horizon at t > bbm~ l as the increase of </> dominates. 

The spectral index n(k) — 1 = dP(k)/d\nk and its /c-dependence dn(k)/d Ink are given 
in Fig. 4. The value of n(k) varies from unity at the larger scales as Ink < 110 to about 
n < 1 at the smaller scale as In A; > 110. This leads to a index-running dn(k)/d\nk < 
in the wavelength range corresponding to the horizon-crossing times during 55 < mt < 70. 
Once again, the physical reason for the running of the spectral index is due to the energy 
transfer from the inflationary background field to the fluctuations as well as the evolution 
of the inflaton field. Although one may refer this energy transfer to the third or even higher 
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FIG. 4: Spectral index n{k) and its running dn(k)/dlnk vs. Ink (k in units of m) with the same 
initial conditions and parameters as in Fig. 1 for the case of N = 4. 

order derivative of the effective potential V"'((j)) of the inflaton, the index-running shown 
in Fig.4 can only be obtained through a proper off-equilibrium dynamics, but not by the 
classical effective potential approach. 



V. THE CORRELATION OF PRIMORDIAL PERTURBATIONS 

A. Correlation function and initial condition of classical field fluctuations 

The equal-time two-point correlation functions of fluctuations have been extensively used 
to study the applications of the O(N) model to the off-equilibrium chiral phase transition in 
relativistic heavy ion collisions where the coupling constant is of order one [15]. The corre- 
lation caused by the off-equilibrium process is generic, because the evolution of fluctuations 
is typically with respect to a time- dependent background. In our case, the mode functions 
of fluctuations f Xt k(t) and f n> k(t) are not trivial plane waves, they depend rather on the 
dynamics of interaction or of self-interaction. Hence, it would be interesting to search for 
the correlations of the primordial perturbations. 

By means of Eq. pip, the two-point correlation functions of the fluctuations are given 

by 

( X (x,t) X (x',t)) = /^^ x - x ')|/ x , k (t)| 2 

1 f A 

= TTi n kdk sinfc l x - x 'll/xk(*)| 2 ' 

Z7T X — X \ J 
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FIG. 5: Mode function |/,r k(*)l vs - k(m) at t = 40, 55, 70 in units of m . All parameters are 
specified the same as in Fig. 2 except the initial conditions Eq. (44) where M?(0) = M^(0) = 
are assumed. 



(71* (x, t)7Ti(x',t)) 



1 



27T 2 |X — X 7 



fcdfc sin fc|x — x'| |/ w k(^)|" 



(46) 



where we have considered that the dynamics of l/^kWI 2 and |/ x ,k(£)| 2 is isotropic in k-space. 

It can be seen from Eq. 1)46)1 that both correlation functions (x(x, t)x(x', £)) and 
(7i"i(x, t)7Ti(x', £)) over large spatial distances |x — x'| are sensitive to the long wavelength 
behavior of |/ Xi k(£)| 2 and |/7r,k(^)| 2 - However, the infrared behavior of the mode func- 
tions during the early stage of the inflation are actually governed by the initial condi- 
tions of the fluctuations as shown in Eq. ()28|). For instance, if we take the initial values 
M^(0) = M^(0) ~ m 2 in Eq. (|Hj). then one can show that with Eq. (|28Jl the correlations 
decay exponentially as exp{ — (m|x — x'|)} = exp{ — (H\x — x'|/2)} when we set if = 2m 
to mimic the constant expansion rate during the inflation. Hence, the correlation between 
the inflationary perturbation modes is negligible if the spatial distance between the modes 
is larger than a horizon size a|x — x'| > H~ x = 0.5m -1 . On the other hand, if the initial 
conditions for the mode functions are chosen as M?(0) = M^(0) = 0, then the solution to 
the mode function can be approximated by 



l/k(f)| S 



3\2 
2> 



,2[i/-(3/2)]JTt 



(47) 



This function is singular as when k — > 0. For such infrared behavior of the mode 
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function, we have 

(X(x,t)x(x',t)) _ (7r(x,t)7r(x',t)) 1 ( , s 

(X 2 ) (vr 2 ) {Ha\x-*\)*' 1 J 

This correlation clearly covers spatial range much larger than the horizon size /T -1 . 

The numerical solutions of \f n ,k{t)\ 2 with the vanishing initial mass terms are plotted in 
Fig. E| Obviously the mode functions given in Fig. have shown very different infrared 
behavior from those in Fig. where we have taken M 2 (0) = M 2 (0) ~ m 2 . However, we 
expect that the tail of the correlation function depends sensitively on the initial conditions of 
the inflaton fluctuations. To illustrate this feature, we plot in Fig. Elthe correlation functions 
(7r(x, t)n(x.',t)) obtained by two different initial conditions: M 2 (0) > and M 2 (0) = 0. We 
see that the respective behaviors of (tt(x, t)7r(x', t)) at scales r = |x — x'| larger than 0.5 vrC x 
(which equals the horizon distance 1/H) are quite distinguishable. With M 2 = initially, 
the correlation function scales as a power law: 

(7Ti(x, t)TTi{x', t)) CX |X - X'|- L8 . (49) 

On the other hand, the correlation function clearly possesses an exponential fall-off with 
distance when the initial condition is taken to be M 2 (0) > 0. This example shows that 
the initial conditions can be recalled from the tail of the correlation function of classical 
perturbations despite that the original quantum fluctuations for those superhorizon modes 
have went through a quantum-to-classical transition during the off-equilibrium evolution. 
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B. Detection of the two-point space-scale correlation 



The primordial density perturbations is governed by the fluctuations of the fields x an d 
7Tj. Thus, the field correlation functions (x(x, t)x( x 'j t)) and (7rj(x, t)7Tj(x', t)) will inevitably 
lead to the correlation of the primordial density perturbation (5(x, t)5(x.', t)). When pertur- 
bations of scale k cross the horizon at a time t, the equal-time two point correlation function 
between (x, t)-(x',t) will yield the correlation between two space-scale modes (x, /c)-(x', k) 
via the mapping formula k = a(t)H [Eq. The density contrast £k(x) represents the 

fluctuations in total energy density at the spatial point x with a scale k. 

Obviously one cannot measure the density (energy) perturbations precisely on a scale 
k and at a spatial point x simultaneously. The essence of 5^ (x) is as follows. When a 
perturbation crosses the horizon at the scale k = aH, the position of the perturbation has 
an uncertainty typically given by the size of the horizon, i.e. Ax = (aH)^ 1 . By virtue 
of the uncertainty relation AxAk ~ 2ir, the scale of the perturbation thus lies within the 
band from k — (1/2) Ak to k + (1/2) Ak where Ak ~ aH. Therefore, <5k(x) describes the 
perturbation in a phase-space cell (x, k) with a size confined by the spatial range from x 
to x + Ax, and by the scale range from k to k + Ak. The volume of this phase-space 
cell is characterized by Ax • Ak ~ 2n. Consequently, (5k(x)5k(x')) governs the correlation 
between two perturbation modes localized within the cells centering at (x, k) and (x', k) in 
the phase (x-k) space. In order to unveil the effects of (x(x, t)x(x', t)) and (vr(x, t)7r(x', t)), 
it is indispensable to decompose the mass density perturbations into the x-k modes in the 
phase space. The discrete wavelet transform (DWT) is designed to do such space-scale [x-k] 
decomposition fl6| . 

In the formulation of DWT, there are two sets of spatially localized bases given by the 
scaling functions |j, l) s , and the wavelet functions |j, 1) TO ; both are characterized by the indices 
j and 1. For a 1-D sample with a spatial size L, the index j = 0, 1, 2, ... stands for a scale from 
kj to kj + Akj in which kj = 27r2 J /L and Akj = 27r2 J /L. The index / = 0, 1, 2 J — 1 denotes 
the location of the spatial point within Llj2? < x\ < L(l + 1)/2 J '. These bases are complete, 
and they satisfy the orthogonal relations s (j, l) s = Siy, and w (j' ,l'\j,l) w = Sjj'Sij'. For 
3-D samples, the DWT bases are given by the direct product of the 1-D bases. Thus, a 
density field \5) can be decomposed into the phase-space modes (j,l) as 

Sj^= s(j,l\6) = /rf 3 x5(x)0 j , I (x) = /^|iW4i(k), 

d*x <KxM,i(x) = jj^y <Wj,i(k) , (50) 

where 0j,i(x) = 0j,i(k) = (fc|j,l) 8 , = and 0j,i(k) = l) w , i.e. 

they are the scaling functions and the wavelet functions in either the ^-representation or the 
/c-representation, respectively. 

Subsequently, the two-point correlation functions in Eq. ()46|) can be rewritten in terms 
of the DWT bases as 

^| /xik (t)|Vj,i(k)^(k), 

(X^v) = /(0il/x,k(OI 2 4i(k)^(k); 
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<WW> = / (^l/^)r^(k)^(k), (51) 

where the time £ is taken to be tj specified by the relation 2ir2 3 /L = k = aH. Since 
a = exp(Ht), one has 



1 

ti = — 
3 H 



j In 2 + In 



(52) 



Thus, Eq. (151)1 can be used to determine the correlations between fluctuations at different 
spatial points 1 and 1', both crossing out of the Hubble horizon at the same time tj during 
the inflationary epoch. 

For example, in the case of a free scalar field \ with mass m, the DWT mode-mode 
correlations can be expressed as 

(2^2(FT^^ l(k *' (k) ' 

p^^T^ €i( k )€i'( k ) • (53) 

Since T/>j,i(k) is localized at k ~ 2tt2 : >/L, and </>ji(k) is non-zero at < 2tt2 j /L, one can 
approximate (fc 2 + m 2 )" 1 / 2 ~ m^ 1 if the j-scales are greater than the Compton wavelength 
of mass m. Thus, 

(XMXl,!') - 7^7^ / d 3 k^(k)^ v (k) = if ^ I' , 

tow) - 7^777, / ^!i(k)fe(k) = o if / ^ r , (54) 

where the orthogonality of the DWT bases with respect to the indices 1, Y has been applied. 
Similar results hold for the field 71*^1. Hence, for the free field case the DWT modes are not 
correlated and therefore, the correlation function (<5k(x)5k(x')) is trivial, as long as the scale 
under consideration is larger than the Compton wavelength of the scalar field. 

However, for the self-interacting O(N) model the DWT mode-mode correlations [Eq. 
(joTj) ] are generally non-zero. Accordingly, the two-point space-scale correlation (<5k(x)5k(x')) 
may become non-trivial and initial condition dependent. Using Eqs. ()39j) - (j43|) . it is straight- 
forward to calculate the normalized two-point space-scale correlation function defined as 

As a numerical example, the normalized two-point space-scale correlations of density 
perturbations under different initial conditions are plotted in Fig. [TJ in which the parameters 
are taken to be H = 2m and L = 2m~ l . From Eq. (|52jh one has 

tj ~ (0.35j + 0.23)m _1 . (55) 

Thus for j =4 — 6 shown in Fig. Q we have tj ~ (1.63 — 2.33)m~ 1 which corresponds to 
the number of e-foldings Htj ~ 3.3 — 4.7. Since the inflation under consideration lasts from 
t — to about t = 70m -1 , the correlations in Fig. [7| actually probe the inflaton dynamics at 
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FIG. 7: The normalized mode-mode correlation functions of density fields in phase space are plotted 
against the separation r (in units of m _1 ) of two perturbation modes with respect to various scales 
j under two different sets of initial conditions. The parameter L is taken to be 2m~ 1 . 

the beginning of the inflation. With M 2 (0) = M 2 (0) ~ m? > 0, the space-scale correlation 
function approaches zero drastically as the distance r between the two modes increases; 
while the correlation deviates from zero and lasts for a long range for perturbations with 
M\ = Ml = initially. 

Figure |H1 plots the same correlations as in Fig. [7|but with respect to j = 8,9, or tj ~ 
(2.75 — 3.38)m _1 , which corresponds to the number of e-foldings Htj ~ 5.52 — 6.8. The 
behavior of correlations under the initial condition M 2 (0) = M 2 (0) ~ m 2 remains the same. 
However, the correlations with M 2 (0) = M 2 (0) = do not coincide with those shown in 
Figure That is, the space-scale correlation of density perturbations in the model with 
M 2 > initially does not experience a significant evolution during the first few e-foldings 
of the inflationary expansion, while the correlation in the model with initial M 2 = does. 
Therefore, the j-dependence of the mode-mode correlation is practical to test the dynamical 
predictions about the evolution (t e -dependence) of the primordial perturbations. Since the 
resolution of DWT basically can be refined as required by the problem in hand, one is capable 
of exploring the physics of the very early universe by means of the mode-mode correlations 
in phase space. 

VI. CONCLUSION 

Using the O(N) model as an example, we have investigated the self-interacting effects on 
the primordial perturbations using a self-consistent off-equilibrium approach. It is known 
that most information of the dynamics and initial conditions of the system will be lost during 
the non-equilibrium evolution. The non-equilibrium evolution can only be probed via the 
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observable remains produced from such a process. In this paper, we have shown that the 
off-equilibrium evolution of cosmological inflation may have two such observable "remains" . 

The first one is the running spectral index of primordial perturbations induced by the 
scale-dependence due to the self-interaction of the inflaton. We found that the running 
spectral index depends essentially on the rate of the particle creation and the energy transfer 
from the background to the inflaton fluctuations as well as the evolution of the background 
field. It is a signature of the energy transferring dynamics during the inflation. Although the 
running index of the O(N) model is small comparing to current data, the negative running 
(Fig. 4) does coincide with the WMAP observation. 

The second remain is the correlation function between phase space modes of the den- 
sity perturbation. Under the influence of the self-interaction, fluctuations created from the 
background field are no longer white noises. Although the one-point distribution function 
of fluctuations in each Hartree mode is Gaussian, the power of the fluctuations is scale- 
dependent, which gives rise to the correlation between the phase space modes of the energy 
density perturbation. Moreover, since the dynamical evolution of the correlation depends 
upon the initial conditions of the inflation, the mode-mode correlation of density perturba- 
tions also provides a window to study the dynamics of the self-interaction as well as the 
initial conditions of the inflation. Here, we would like to emphasis two important results. 
First, the initial-condition-dependency of correlation functions is irrelevant to the number of 
fields in action. Second, the inflationary density perturbations after the superhorizon evolu- 
tion are not the fluctuations of thermal equilibrium state. Therefore, the dependence of the 
correlation function upon the initial conditions does not contradict with the reheating of the 
universe, which is generally produced by the interactions between and other fields. Thus, 
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we may expect that the non-trivial mode-mode correlation in the phase space is detectable 
via a DWT analysis on the CMB temperature map, or other observations on the large scale 
structure relevant to the density perturbations [l8| . 

Although, the O(N) model is just an example to illuminate the off-equilibrium effects 
from the self-interaction of the inflaton field, we believe the implication from what we have 
found in this paper is also useful to other interacting inflation models. For instance, the 
scale-dependence of perturbations drawn from the self-consistent off-equilibrium dynamics 
is generally different from that obtained by the third or the higher derivative of a classical 
effective potential. Therefore, the running spectral index determined by the formalism of 
various effective potentials would be questionable if the interaction or self-interaction of the 
inflaton is substantial. 
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